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Homogeneous isotropic spatial flat cosmological models with two torsion functions in vacuum 
are built and investigated in the framework of de Sitter gauge theory of gravity. It is shown that 
by certain choices of parameters of gravitational Lagrangian the cosmological equations have some 
, exact constant solutions that turn out to be stable de Sitter critical points of dynamical systems and 
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can explain observable acceleration of cosmological expansion. The role of the space-time torsion 
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provoking the acceleration of cosmological expansion is shown. 
PACS numbers: 04.50.Kd, 98.80.Jk, 98.80.-k, 11.15.-q 



CJ ■ I. INTRODUCTION 

cr: 
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bX), Recently there has been a burst of activity dealing with quadratic gravitation. For example, the curvature- 



squared terms added to the usual Einstein action with cosmological constant have played a role in two 
recent investigations of four-dimensional gravity: in critical gravity [1], and in a pure Weyl-squared action 
considered by Maldacena [2]. 

The critical gravity provides a consistent toy model for quantum gravity as a useful simplified arena for 
studying some aspects of a potentially renormalisable theory of massless spin-2 fields in four dimensions. 

The conformal gravity theory has been advanced as a candidate alternative to standard Einstein gravity. 
As a quantum theory the conformal theory is both renormalizable and unitary, with unitarity being obtained 
■ because the theory is a PT symmetric rather than a Hermitian theory. Because the variation of the conformal 



action leads to fourth-order equations of motion, it had long been thought that the theory would not be 
unitary. However, as has been shown by Bender and Mannheim [3] that one can find a realization of the 
theory that is unitary. Consequently, conformal gravity is to be regarded as a bona fide quantum gravitational 
theory. The conformal gravity theory can quite naturally handle some of the most troublesome problems in 
physics, the quantum gravity problem, the vacuum energy problem, and the dark matter problem. [4] 

As a modified gravity theory quadratic gravitation has been used in cosmology [5]. In order to explain 
observable acceleration of cosmological expansion some authors introduce torsion terms in quadratic grav- 
itation [6]. The quantum aspects of torsion theory and the possibility of the space-time torsion to exist and 
to be detected have been discussed in [7]. The astronomical observations show that our universe is prob- 
ably an asymptotically de Sitter (dS) one with a positive cosmological constant A [8]. If a gravitational 
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theory of Yang-Mills type is constructed starting from de Sitter gauge invariance principle, its gravitational 
Lagrangian naturally turns out to be the one of quadratic gravitation with torsion as will be shown in this pa- 
per. Therefore, a investigation of quadratic gravitation with torsion and its cosmological solutions expressed 
by de Sitter critic points will be carried out. The field equations will be derived. These equations are quite 
different from the equations obtained from Riemannian geometry based quadratic Lagrangians when varied 
with respect to the metric. Applying to the space flat FRW cosmology some de Sitter critical point solutions 
will be obtained. The stability of them will be analyzed. 

The paper is organized as follows. In section U, starting from a CUfi'ord algebra C (3, 1) the gravitational 
Lagrangian of a de Sitter gauge theory is constructed, the Lagrange equations of gravitational fields are 
derived. Applying them to a spatial flat universe the cosmological equations are obtained in section m. The 
vacuum solutions of these equations in two specific cases are presented in section TV. These two models 
correspond to the conformal cosmology of Mannheim [4,9] and the zero-energy gravity of Deser and Tekin 
[10], respectively. In contrast to them, the tetrad and the spin connection are taken to be the basic field 
variables and the torsion plies a important role here. In these specific models the cosmological equations 
are written as some dynamical systems, the real de Sitter critical points of them are obtained. Among these 
points, the stable ones which turn out to be exact constant solutions and describe the asymptotic behavior 
of the universe are found. In section V some concluding remarks are given. In Appendixes the calculations 
for stability analysis are presented. 



II, LAGRANGIAN AND FIELD EQUATIONS 

We begin with a brief introduction of a de Sitter gauge theory. In a gravitational gauge theory coupled 
to matter sources involving Dirac fields it is convenient to take Dirac matrices y/ and their commutators 
CT/y - 2 [77,77] as the basis of the gauge algebra. In this case we are led to a de Sitter gauge theory. 
Let {7/} (/ = 0, 1, 2, 3) be a basis of an inner product space with signature (-, +, +, +). A Chfford algebra 
C (3, 1) can be constructed by introducing the condition 

jijj + jjji = ^mji- (1) 

with 77// =diag(-l; 1; 1; 1). There is a 10-dimensional subspace of C (3, 1) which is a Lie algebra with basis 
757/ and crjj = ^ [7/577]- This is the Lie algebra of a de Sitter group. We can introduce a connection 
[11,12] 

oj^T+ y75e, (2) 
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defined by 

t-e\yi®dx^, (3) 

and 

where / denotes a constant with the dimension of length. The curvature of oj is 

n-^/a; + ^[w,a>] = R+yr5T-iv, (4) 

where 

R = jr + i[r,r], 

T = Je + [r,e], 

V = ^[e,e]. (5) 

The Lorentz curvature R, the torsion T, and the cosmological term V are given by, respectively, 

R = lR'-'f,y(rij®dj<f' Adx", 
8 

T = ^T^^yo-ij^djd' Adx\ 

V = e'ij,e\(rij ®dx!^ h dx", (6) 



with 



R"^.v = df,r'\ - dyT"^ + 7]KLr"'jjr'^\ - riKir^'yT'^^, (7) 

and 

T'^^y = - dye'^ + r'j^e'y - T' jye\. (8) 

Based on the local gauge invariance principle the gravitational Lagrangian can be made up of a quadratic 
term of the curvature D and its Hodge dual *Q: 



(9) 



X = -irr (*Q A Q) = |l/?^/-/?/'v - ^r^r^pr/'' + ir^/? - iir j e, 

where 

e = det|e^^|. (10) 
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In four dimensional spacetime the Gauss-Bonnet term -^Pg [RuvZtR'^^'^^ - ^-RfiyR'^^ + is purely topo- 
logical and then the Lagrangian can be taken as 

£ = -^Tr(*n A Q) = V - - + ^r^R- i2rj^. (ii) 

For the sake of a neater argument we extend the Lagrangian to including the coefficients 



and rewrite (10) as 



£ = i^r^R^yR^" + aF^R^ + yV^T^'ypT + ^r^R - lir^^ e = Le, 



(13) 



with 



L = pr^Ri.yRi'^ + al'^R^ + jl'^T^ypT/P + ^I'^R - 121'^. (14) 

£ is just the Lagrangian of quadratic-curvature gravities [10] with torsion. 

The variational principle yields the field equations for the tetrad and the spin connection 



^'^ = esi/, (15) 



5F"^ 



where E'^j and are energy- momentum and spin tensors of the matter source, respectively, the variational 
derivatives are given by 

6£ 

= (le^^RP^Rp^ + le^pRP'^R'j^^ - e' ^Rp^RP") + al'^ {46" Ry^, - e\R)R 
+yr^ (Ae'''T'y,T;^,' - Ady (e'^Tp/) - e^T^^T r + {Ae'^Tp/) e'^rdyCK') 

^e'^'Ryp - ^e>J + I2r^e'p]e, (16) 

= {iprWleW/ - ej'dyR/ + (ejV - e/'R/) e^rdyeK^ 
+e/TV/?/ + ei'Ty.R/ - eii^ry.R/ - eI^T^'yrR/] 
+2al-\{ei'e/ - ejW) T^rR + {e/e/ - e/e/) (Y^^yR - dyR) 
+ {ei^'e/ - eiW)Re''rdyeK''i + Ayr^ejye/r^ 

+h-\ejW - ejW) r^vr + (ej'e/ - efe/) {v\y + e^,dyeK')\}e. (17) 
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That may be, the two main field equations are rather complicated. They really look nothing hke the famihar, 
well-analyzed equations of GR. To help understand the significance of these new relations, and to use 
our previous experience, we will do a translation of (16,17) into a certain efi'ective Riemannian form- 
transcribing from quantities expressed in terms of the tetrad ef and spin connection T^^ ^ into the ones 

expressed in terms of the metric g^y and torsion T'^^y (or contortion ^"'^y). 



As is well-known, the affine connection Y'^^y can be represented in the form 



rV = ei^d/y + e/e\Y\f, 

= {/v] + K\y, (18) 



where \ix^v\^ K'^nv are the Christoffel symbol and the contortion, separately, with 

I^A. ^ / 'T'yl -I- T" _|_ T* ^\ 

^ //y — \^ fXV ~^ fXV VjU j ■> 
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T nv — e/'T ~ r ~ r yp. (19) 



Accordingly the curvature can be represented as 



+ {/J K\y - {/y} ^rV;. + K% - {^-^j KP^y, (20) 



where 



is the curvature of the Christoffel symbol. 

m. COSMOLOGICAL EQUATIONS 

For the space flat Friedmann-Robertson- Walker metric 

g^y = diag (- 1 , a {tf , a {tf , a (tf) , (21) 



we have 



{o%} = 0,{o°4 = {,%} = 0,{,«,} = aa5,-,-, 

{o'o} - 0, {/o) = {o',) = -S), {/k} = 0, /, j, k, ... = 1, 2, 3. (22) 
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The non- vanishing torsion components with holonomic indices are given by two functions h and / [13]: 

7^1 10 - T220 - T^3Q - a^h, 

Tm - T231 - T^n - a^f, (23) 

and then contortion components are 

K^n — K^22 — K^22 — o^h, 

^^23 - ^^31 - ^^12 = ~2^f' 

K\2 - K^n = K\i = ^af. (24) 

Among the torsion components, only the pseudotrace axial ingredient given by / couples to spinors in a 
minimal way. The scalar mode h of torsion could be considered as a "phantom" field, at least in the matter- 
dominated epoch, since it will not interact directly with matter; it only interacts indirectly via gravitation. 
The non- vanishing components of the curvature R'^^a-fiv are 

R^ioi = /?S202 = /?f)303 = [h + Hh+ hj , 

R^]123 = -R'^^213 = ^^]312 = a^f{H + h) , 
^{}203 - -^{)302 = ^f)301 = -2^iHf+ fj, 

R\}2i2 = R\}3n = ^323 - fl' {h^ + 2Hh + h^- ^/J , (25) 
/?„oo = -3 H -3 h -3H^ - 3Hh, 

R{]n = (h +3H^ + 5Hh+ h +2h^ - ^/j , (26) 

R{] = 6H +12H^ + imh + 6h +6h^ - ^f, (27) 
where H -a (t) /a (t) is the Hubble parameter. Using these results and (16 — 20) we can compute 



eio-^ = {(^ + 3a) [-n I^H + hf -24 i^H + h^H(H + h) 



+\2h (h + 2H) (h + Hf-6{h + Hf f + 

+y[nh^ + 6f) + 3H^ + 6Hh + 3h^ - jf - nr^}e, (28) 



en-^ = -rV{0S + 3a)[-4(77 + - %(h + h^(H^ + Hh) 

+4h (h + 2H) (h + Hf-2{h + Hf f + 
-27(2 h +SHh + ft^ + + 2 + /i) + 3H^ 



_6£ 

<5roi 



+4Hh + h^- - ni-^}e, (29) 



+ (5/3 + 18a) (Hh+hh +h^H^ + (/? + Sof) \2h^ -ff -\hf\ + ^h}e. 



(30) 



5£ 



= r V{2 (jS + 6a) (77 + ii) + 6 (;8 + 4a) //^ 



6T'h 

+2 (5;8 + 18a) ///i + + 3a) [Ah^ - f) 

(31) 

Suppose the matter source is a fluid characterized by the density p the pressure p and the spin ^/y'^. The 
system of field equations (15) consists of four independent ones: 



6Z S£^ 
6£ 6£^ 

6£ 5£^ J J 

= ei SQi , 



Using (28-31) the Lagrange equations (32) can be written as 

v2 



^01 J gp)\ 

S£ ^■^'P 3 /ooN 

- ^3 Si2 ■ (32) 



{J3 + 3a)[-n(^H + hj -2a[h + h^H{H + h) 
+ \2h (h + 2H) (h + Hf-6{h + Hf f + 

+7(18/?^ + 6f) + 3i7^ + 677;i + - jf - nr^ - fip = 0, (33) 



2 

(yS + 3a) [-4(h + - s(^H + hj(H^ + Hh) 
+4h (h + 2H) (h + Hf-2{h + Hf f + ^-f\ 
+27(2 h +mh + +f^-2(H + hj- 

-AHh -h^ + + 12/-^ + fp = 0, (34) 
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-2{0S + 6a) + /j) + 3 08 + 4a) (hH^ + 2HH +2h Hj + (5/3 + 18a) (h h +h h +h^H^ 

+ OS + 3a) {ih^ - / / - ^/i/') + ^/i} - ;'^o/ = 0, (35) 

f{2 (fi + 6a)(H + hj + 6(p + 4a) 

+2 (5/3 + 18a) Hh + (J3 + 3a) (4h^ - f) - 47 + ^} - fsn^ = 0. (36) 
Assuming s^^v'^ = (i.e., the source spin current is negligible), the Eq. (36) reads 



/{2 OS + 6a) (// + /i) + 6 OS + 4a) 

+2 (5/3 +\Sa)Hh + (J3 + 3a) (4h^ - f) - 47 + ^} = 0, 



(37) 



and gives 



or 



/ = 0, (38) 
20s + 6a) + ij) + 608 + 4a)H^ 

+2 (5yS + 18a) Hh + (J3 + 3a) {4h^ - /) - 47 + = 0. (39) 



Therefore, we have two cases. 

In the first case, / = 0, the Eqs (33) and (34) read 

\2 



(fi + 3a)[-(H + hj -2{^H + h^H{H + h) + h{h + 2H)(h + H)^] 

+-yh^ + -H^ + -Hh + -h^ - - — = 0, (40) 
2 4 2 4 12 



and 



08 + 3a) [-in + hj - 2 (// + hj (h^ + Hh) + h(h + 2H) (h + H)^] 
+^7(2 h +8Hh + h^y^iH + hj-^H^-Hh-'^ + 31-^ + ^ " ^' 



which lead to 



(41) 



72 

(2y -l)h -2H^ + (87 -3)Hh- {2y +\)h^ + 8/"^ + - Oo + 3;?) , (42) 

6 
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and 



-4j^h +y^4H^ + Sil-4y)Hh + 4i2y+\)h^-^l^(p-3p)-^^h 
+ \6H^hy + (287 - 64y^) H^h^ + 8 (47 + 1) yh^H - 4 (7 + 1) yh"^ 

1 /16 1 , \ 9 (6r + 1) 7 



120S + 3a) 3^ ^ 36 ^ {p + 3a)f 

= 0. (43) 

So we have the equations (42), (43) and 

(j3 + 6a)(H + h^ + 'i(J5 + 4a) {hH^ + 2H H +2h 7/j + (5^ + 18a) h +h h +h^H^ 

+2{p + 3a)h^ + ^h^Q, (44) 

for the unknown functions H and h. 

In the second case, / satisfies the condition (39). The Eqs. (33) and (34) yield 

1 P- 

(2y -l)h - 2H^ + (87 -3)Hh- {2y + l)/?^ + -f + 8/"^ + _ (p + 3p) , (45) 

4 o 



and 



-472 h +7 [4H^ + 8 (1 - 47) //^i + 4 (1 + 27) h^-f- ^fiB -y^h 

+ l6H^hy + 47 (7 - I67) H^h^ + 87 (1 + 47) h^H - 47 (1 + 7) - AyHhf + yh^f 

/16 1 , \ , 1 , /32 2 , \ 1 

, (_ , , 3;,)) i.^ ^^//H (1 - 47) . -/^«) . ^^//Z. 

.f(1.27)/.^4(1.27)/^(p.3p)/.^.^|l^/.^-(l.l/^(p.3,))/ 

^ 16^^ + la/ - - 12 3af " ^ + ^p) - 1/^ (p + 3pf - ^ 

= 0. (46) 



The Eqs. (45) and (39) gives 



r = ^y^'^-^h +4H' + S\l+4y^^^j-^jHh + ^4-8y^i^-^jh 

^1-87 ^ 32^ ^ ^^lMzV,3p). (47) 
13 3/3 ^ ^ 



10 



Substituting into (45) and (46) yields 



-I2y 



2/? + 4a -2 



+ |-32| {y(J3 + 6a)-2(J3 + 3a))Hh + 8^ {y(J3 + 6a) + 2(J3 + 3a))h'^^ h 
, / 8y-l / y06 + 6a) , ^\ 4 17^+48a 17^+48a . ^ 

P P P 



13 + 3a 



Ij 



24y/3 + 96ya - /3 - \2a 



- 384y 



P + Aa B + 4a 



+ 



/3(j3 + 3a) 

-5/3 + \2yp + 4Sya - 6a , ^6 + 4a ^ ^ /3 + 4a ,2. , ^ ^ 
-y + 96y h 2y r (p + 3p) 



-f (p + 3p) 



Hh 



^l-8y 8y-l ^ 48yj8 + 192yQr - T^S - 24a ^^2-'^'^" 



p \6{l3 + 3a) 



p{fi + 3a)f 



12(j8 + 3Qr) 



■p + 



8y - 1 yS + 4a ,2 
8 pip + 3a) 



= 0. 



(48) 



and 



l-8y ,,j6 + 3a ;8 + 3a,2 



4;S 



+ 16- 



pf 



+ 



3;S 



-l'{p-3p) 



(49) 



Differentiating (47) gives 

• (iS + 6a) •• / {p + 6a)\ ■ ■ ( (P + 6a)\ ■ I (P + 6a)\ ■ 

-ff = -4y ^ ^ ^ h -4 l+4y '^ ^ \h H -AH H -4 l+4y ^ ^ ^ // /t - 4 - 8y ^ ^ ^ U 



;8 

(j8 + 6a) 



■(p+3p). 



Substituting it and (47) into (35) and letting sq\^ =0 give 

iP + 3a) 



77 + 1 - 4y- 



8y 



h+2HH+2\l-^{p + 3a)\hH 



P 



P 



p(P + 6a) 4p 
{p +3 



pf 



3p 



(P + 3a) ^2 



3p 



= 0. 



(50) 
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So we have the equations (48), (49), and (50) for the unknown functions H and h. The unknown function 
/ is given by (47). 

IV. TWO SPECIFIC MODELS 

In order to emphasize the geometrical nature of the effect of acceleration of cosmological expansion we 
concentrate on vacuum solutions in two specific cases and discuss only the acceleration solutions. 

A. When ^ = -3a 

This corresponds to conformal (Weyl) gravity which has been investigated by numerous authors [recent, 
see 2 and 9] but it must be pointed out that the principle and structure between the theory here and higher- 
derivative gravity in Mannheim's theory are quite different. 

According to last section, the equation (37) gives two cases. 

In the first case / = 0, the functions H and h now satisfy the equations (40), (41) and (44), i.e., 

i6y+l)h^ + + 2Hh - 41'^ = 0, (51) 

(47 -2)h-2H+ (I67 -4)Hh + {ly -l)h^- 3H^ + 12/"^ = 0, (52) 

H + h+2{H + h)H 

+ {H + h)h +hH^ + h^H + -^h = 0. (53) 

12Qr 



Eq. (51) has the roots 



-H± V-6y//2 + 4(6y+l)/-2 
(67 + 1) 



Eq. (52) gives 



and then 



1 (87-2) 1 , 3 , 61-^ 

(27 -1) (27 -1) 2 2(27-1) (27 -1) 
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Substituting (54), (55) and (56) into (53) yields 

•• _ ( 4Sy^ - 50y^ -ly + 2 _ 87 - 1 yJ-6yH^ + 4(6y + I) l'^) ■ 2 (5047^ + 3247^ - 267 - 3) 
" " ~[ 27(27- 1) (67 +1) (67 +1) (67+1)3(27-1) ' 

8407^ - 472 - 67 - 5 27-1 



(2y - 1) (67 + 1)^ 7/2 24a7 (67 + 1) 
.,2 _l1*:qo„4 _L-7/i/i„3 



47(67+1)^(27-1) 7(67+1)2/2 24a7(67+l)/> 



Let 

We have the dynamical system 



H ^ X, 

( 4Sy^ - 5O72 -77 + 2 _ 87 - 1 V-67//2 +4(67+ l)r2^ 
X — —\ — — — T—H 



27(27- 1) (67+ 1) 



4y 



(67 +1) 



X 



+AH^ +BH + {CH^ + D) ^-677^2 +4(67+ l)/-2 



(58) 



where 



The Jacobian elements are 



A = 



B = 



C 
D 



2(50473+32472-267-3) 

(67+1)3(27-1) ' 
84O73 - 472 - 67 - 5 



27-1 



(27 - 1) (67 + 1)2 7/2 24a7 (67 + 1) ' 

-47672 + 25927^^ + 74473 - I87 + 1 

47(67+ 1)3(27- 1) 
3 (27 +1) 27-1 



+ 



7(67+ 1)2/2 2Aay(6y+l) 



^ = 0^ = 1 
dH ' dX ' 



(59) 



dX _ ( 4873-5072-77 + 2 , 3(87-1)// ' 

dH [ 27(27- 1) (67+ 1) ^ 2(67+ 1) V^^67^2T4(67TTyF2, 

+A//3 +BH + (CH^ + D) ^-67/^2 + 4(67+ l)/-2, 



X 



dX 
dX 



4873 - 5O72 -77 + 2 



H 



87-1 



27 (27 - 1) (67 + 1) 47 (67 + 1) 



^-67/^2 + 4(67+ l)/-2. 



(60) 



13 



The critical point equations are 



X = 0, 

AH^ + BH + (CH^ + d) ^-67772 +4(6y+ 1)/-2 = 0. (61) 

Rationalization gives 



+ aH^ + bH^ + c^O, (62) 

where 



a - 



b = 



2 {f-AB - nC^y - 2C^ + 6CDjfi) 

(A2 + 6C2y) P ' 
s2/2 _ 4SyDC - SCD + 6D^jf 
(A2 + 6C27)Z2 ' 



, . -j^^^D^ (63) 

(A2 + 6C27) /2 



The equation (62) has the roots 



where 



and 



„f.(_|.VS)"%(-|-VS)'"-| 



«|.(_|.VX)"'..(-|-VAfl|, 

«| = (-|,,G)"'^2^(_£_Va)'".-| (64) 



2 . 1 
q = —a - -ba + c, 
^ 27 3 



A^^^)%(ff, (65) 



-1 + V3/ 

ta = . (66) 



14 



Now we have the critical points 

In order to analyze their stability we give the parameter a and y specific values and then obtain the 
results: 
When 



the equations (62) become 



1 ,2 1 



A'iXH^f - UlOQH'^f + \519m^l^ + 3600 = 0. 



It has the roots 



W « 1.4024/^,77^ « -. 19478/r - 2.0 X 10'^ i/f,H^ « 10.596/r + .92212?-/^, 



the first root = I. AOlAj corresponds a positive critical point 



At this point, for 



1.1842//,X = 0. 



-H + V-6y7/2 + 4 (6y + ifj^ _ 1.725 

(6r + 1) " / ' 



the dynamic system (58) reads 



The Jacobian 



H 

X 



X 



3 p ^3 P yi f 



430.76 2. 325 
/2 / 
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has the eigenvalues -1.1625// - 20.7222?//, -1.1625// + 20.7222?//. Therefore, the critical point 



is stable, where 



1. 1842//,Z = 0, 



, 1-725 , ^ 



When 



the equations (62) become 



miH^f + 5226H'^f + 2519H^f - 2312 = 0. 



It has roots 



« 0. 4412//^ « -1. 6508//^ + 0. 37826^//^ 77^ « -1. 6508//^ - 0. 37826^7/^ 



the first root - 0. 4412//^ corresponds a positive critical point 



H = 0. 66423//, Z = 0. 



At this point, for 



-H- yj-6yH^ +A{6y+\)l-^ 1.488 



(6r + i) 



/ ' 



the dynamic system (58) reads 



H = X, 



X ^ - 



17 1 / o 40 



X 



596 , 692 .184. 136, 

-rBH + ( + — r-) 

125 75/2 ^125 75/2^ 



V 2 /2 



The Jacobian 



M = 



1 

10.365 3.4807 



has the eigenvalues -1.74035// - 2.70854///, -1.74035// + 2.70854///. Therefore, the critical point 



77 = 0.66423 //,X = 0, 



16 



is stable, there 



1.488 
^,/ = 0. 



For the two points 

X ^H^ 0, 

which corresponds to a de Sitter spacetime. 

Following Lu and Pope [1], we chose = which means 

"^-48- 

In contrast with them we deal with a de Sitter spacetime with torsion and the gravitational Lagrangian 
including a term yr^T^ypT//'^. When we chose 



(59) and (63) give, respectively. 



and 



, 508 „ 156 ^ 412 „ ^ 

A = ,B = — —,C = ,D = 0, 

3 f 3 



17000 , 27378 

a 7r,b = T,c - 0. 

431/2 431/4 



The dynamical system (58) becomes 



H ^ X, 



X = 



508 . 156 412 2 /3 ,2 



(62) becomes 



431/2 431/4 ' 



and has the roots 



Hi=0,H2 = 



8500 - 103 V5698 „ „ 8500 + 103 V5698 



431 



431 



//. 
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Therefore we get three critical points 



At the point 



the Jacobian 



has the eigenvalues 



This point is a center. 
At the point 



the Jacobian 



has the eigenvalues 



Hi = 0,Zi = 0, 
H2 = 



8500 - 103 V5698 
431 



8500 + 103 V5698 
431 



Hi = 0,Xi = 0, 



ll,X2 = 0, 



//,^3 =0. 



^0 1 ^ 

156 n-6_^ 
\ f ^ I ) 



1 

180.44 4.7728 



-2.3864// + 13.2191///, -2.3864// - 13.2191///. 



This is a stable critical point, where 



= 1.1472//,/ = 0. 



At the point 



the Jacobian 



M 



1 

84.5 46.4 
^ f I 



has the eigenvalues 



This is a unstable critical point. 
If we chose 



-48.15//, 1.755/Z. 



(59) and (63) give, respectively. 



596 _ 164 _ 184 _ 112 
" ~125"^ " ~25f'^~ 125'^ ~ 25P' 



and 



474 , 459 224 



203/2'" 203/4'^ 29 fi' 



The dynamical system (58) becomes 



H ^ X, 
X = 



17 2 



3~ 101 



596 , 164 184 , 



(62) becomes 



It has a real root 



At the critical point 



474 459 224 

^ + t//^ jH^ T=0 

203/2 203/4 29/6 



1.30134//. 



1.30134//,X = 0, 



the Jacobian 



^0 1 ^ 

36. 265 3.3321 
V P I 
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has the eigenvalues 

-1.666/Z - 5.787///, -1.666// + 5.787///. 
This is a stable critical point, where 

= -1.613//,/ = 0. 

In the second case, the functions H, h and / satisfy the equations (33-36) wiiich now read 



2 4 , 8 4(6r+ 1) , 16 

f = + Hh + -h^ T, (67) 

l-Sy 1-87 1-87 (l-87)/2 



-2H +2(2y-l)h -3H^ + 4 (47 - 1) Hh + (2y -l)h^ + ^^^f + 12/"^ = 0, (68) 

H + h +2(H + h) H +(H + h) h +hH^ + h^H+-^h^O, (69) 

12a 

H + h+H^ + Hh-^ + ^0. (70) 
5a 12a 

Eqs. (67), (68) and (70) give 

u 4(87-3) ^, , 2(4y + 3) ,2 2(167-1) , 87-1 „ 

h— n — Hh H — n + — — U. 

87-1 87-1 87-1 7(87- l)/2 247a 

87 + 3„2 . 24y-ll „, ^47 + 3,^ , . I67-I . (87-l)(27-l) 

^-^—r"''-^-^'' +%(87-i)/2+ ^ ' 

and then 

8 ^ • -3 + 87, ■ -3 + 87„- 4(3+47) • 

87 -1 87 -1 87 -1 (-I+87) 

Substituting into (69) yields 
This equation and (70) lead to 



+ hH^ + h^H + -^h = 0. 



h^O. 
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Then (69) becomes 



It has the solution 



The deceleration parameter is 



When 



we have 



and then 



H +2H 0. 



-H^ + C, 



H= ^fc 



g2 VC(/-;,i) ^ Y 
g2Vca-Jo) _ 1 ' 



H . C 



c = o 



dH , 
K-dt, 



1 1 



B. When /3 = -4a 



In this case the gravitational Lagrangian is the square of the traceless Ricci tensor R/^y - R^^y - \gfiyR 
[10]. 

According to section HI, the equation (37) gives two cases. 

In the first case / = 0, the nonvanishing functions H and h satisfy the equations (40), (41) and (44), i.e.. 



H + h] +2{H + h]H(H + h)-h(h + 2H)(h + Hy 



+ 



67+1 
4a 



1 



1 



1 



+ —W + —Hh - -r^ = 0, 



4a 



2a 



a 



(71) 



4^H + hj-Syh +4(27 + 1)^^ -4(87-3)^177 + SH^ - 321'^ = 0, 



(72) 



They can be rewritten as 



and 



Let 



We have 



and 



1 , 1 , 1 
2y 2y loay 



2y + 3, 5y-2\- /By- 3, 2„, 

H = — ^- h+— H\h+\-h: h--H\H 

\ 2y y I \ 2y y 

2y 2y loay 



4/ h -4y^i2y +\)h^- (87 -2)Hh + H^-^^h 

+47 (1 + 7) ^j"* - 87 (1 + 47) h^H + 4y (I67 - 7) h^H^ - l6yH^h 

+ /_ i6(2r^ , 5z±iU^ . (32^2:^ . i-1 HI. 



P 4a I \ f 2a 

P \ f 4a /4 a 



= 0. 



I 27 27 / \ 27 7 

1 7 1 , 1 
+—h^H+—h'--—h, 
2y 2y loay 

/ 27 + 3 57-2 \ /87-3 2 , 

\ 27 7 / \ 27 7 

27 - 1 27 - 1 T 2y — 1 

27 27 loa7 



, /(27+ 1) , 2(47-1) 1 , 8 \ 

\ 7 7 7 7PI 

r r 7 7 
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, 4 1 \ o 16 1 , 



= 0. (79) 



The constraint equation (79) has the roots 



where 



I'^-^i aI(^) -c^Y(H,h), (80) 



\ r 7 7 r^v 

(y+l) ,.4 2(4y+l) ,3„ , (16y-7) ,2„2 4 3 

c = n n H H n H -rhH y 

y y y y 



y'^V- \6ay^ j \ 7^P 8ay 

^ 16^)^' ^ y^ ~ W^"'' ^^^^ 

So we are left with only three independent unknown functions h, H, and X, which satisfies the equations 

H = X, 
h = YiH,h), 

X . (JjLtlH.'-l^H]YiH,h)J'-I^h-'-H]x 



2y 7 I \ 2y y 

2y 2y loory 



The critical point equations are 



The Eq.(86) means 



(84) 
(85) 

^Zzlh^H^^h'-^h = 0. (86) 
2y 2y 



h^O, (87) 



X 


= 0, 


Y(H,h) 


= 0, 


27-\ 


= 0. 


\6ay 



or 



hH + h^- — ^0. 



(88) 



For 



the equations (85) has the roots 
So we have the first pair of critical points 

V -0 h-Ci H - + \ ^ 



For 



the critical point equations become 



X ^ 0, 

1 



5a lOOya 5yP 

1 1 16r-l 16 
+ 1 ^ + ^ + ^-t;- + 



100a2 320ya^ lOOy^fa 25yH 
8r- 1 1 



25600q;372 400y^fa^ 

= 0. 

The equation (91) has the roots 

.;.(_|.VSy'%(-|-VKr-f, 
.l = (-|.^G)'^.(-f-VA)'".^ 

where 



A = 


(fr^(f)^ 


OJ - 


i(-l+V3i), 


P = 






2 , 1 

—A^ - -AB + C, 
27 3 
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with 



,1 3 8 \ A 1 1 8 

A = - — - — — + —T ,-- = +— - — — + 



5a 20070; 57/2/' 3 I5a 200ya l5yP 

1 1 I67-I 16 
B = + + + 



100^2 3207a2 IOO72/2Q; 2572/^ 

C^-l + ^ ] (94) 

\25600q'372 AOOy^Pa^l' 

The equations (89), (90) and (92) give the critical points {X,H,h}. Every one of these point corresponds to 
a de Sitter spacetime. 

The dynamical system (83) has the Jacobian elements 

^ .0^=0^ = 1 
dH ' dh ' dX ' 

dH ~ dH' dh ~ dh' dX ~ ' 

dX 57-2 , , / 27 + 3 57-2 \dY 2 27-1 2 
■Y(H,h) + { — - — h + — H\- X+ '-^ 



dH 7 \ 27 7 I dH 7 27 



where 



dh ly \ 27 7 I dh 2y 

2y - I 27 - 1 , 2y — \ 

+— hH + 34; h^ - -i- , 

7 27 16a7 

|I = !21Z2,_2„, (95) 

dX 2y y 



dH 7 7 

^ ^/iH ^h^H + ^hH^ + \h^ - -i^h - -^H I , (96) 



2 yy^ y2 y2 y2 872a 80:72 



dY {2y+l), (47 - 1) „ 

:sr ^ ^ 

dh y 7 



1 / 1 . 3 , 3 , 1 . 67 + 1 1 



In order to analyze their stability we give the parameter a and 7 specific values and then obtain the 
results: 

For the critical point X - 0,h - 0,H - 2/1, corresponding calculation indicates it is unstable for 

or - r^f. 

In the case X -0, hH + h^ - - 0,v/e have 



, ^h' + -^h'H + ^hW + - -t^h - r-^H . (97) 

2 \7 7 7 7 oy a oy a 
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When 



the equations (92) and (90) give 

h 
h 



1 - 

2 - 

hi = 



1.9814 1.4076 1.4341 

4.4493 + 3.3485? 



4.4493-3.3485? 



At 



for 



1.4341 
'l 



2 -<^ = °- 



the dynamical system (83) has the form 



H ^ X, 
h = Y(H,h), 

X = -(7h + 3H)Y(H,h)-i7h + 3H)X 



with 



Its Jacobian 



Y ^ 3h^ + 2H^ - ^ 



128 



1 

21.324 12. 665 n 
303. 83 185.26 14.288 



has the eigenvalues: -0.39225// + 11. 048?//, -0.39225// - 11. 048?//, -0.8385//. The critical point 

/1.4076\ „ /1. 4341 
h\ - I ; — \,Hi - 



I 



I 
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is stable, where f -0. 
For 



11.886 




the dynamical system (83) has the Jacobian 



M = 



I 1 

9. 8509 4. 2259 q 



/ 

173. 12 
P 



I 

17.401 
P 



14.288 



with the eigenvalues: -22. 33//, 6. 1339// + 1. 6776?//, 6. 1339// - 1. 6776?//. The critical point is unstable. 
In the case / 0, (33), (34), (35) and (39) read (in vacuum) 



[h + hj + 2(^H + hjH (H + h) 

-h{h + 2H){h + Hf + ^{h + Hff - ^/ 

6y+l-, 1 1 8y-l-, 1-, 

+ —H^ + —Hh + -^T^f - -r^ = 0, (98) 

4a 4a 2a loa a 



(H + h)-(Hh+hh +h'H^ - ^ {2h^ -ff -hif^ + " 0' (99) 

4a (h + J^- 4aHh - a {4h^ - /) - 47 + ^ = 0. (100) 

/2 = 4 ^7/ + /jj - 8y /? +4 (27 + 1) /i^ - 4 (87 -3)hH + SH^ - 321'^, (101) 
These equations have the solution 



1 1 



^ = (64-4^'' ^'''^ 
^ ^'^y (103) 



(I67- 1)Z2' 
32(87- 1) 
'/2(167- 1)' 



f -liglzil. (,04, 



For 



7 > — , 
^ 16 
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or 

r<o, 

we have a de Sitter solution 



When 

Irl » 1, 

we have 




a value speculated [14]. When 

1 1 
16 8 

/ is real. In this case, it is the pseudotrace axial ingredient / of torsion that produces the effect of accelera- 
tion of cosmological expansion. 

V. CONCLUSIONS 

Stating from a de Sitter gauge theory a gravitational Lagrangian (13) which is identified with the La- 
grangian of quadratic-curvature gravities with torsion has been constructed. The cosmological equations 
(33-36) for spatial flat universe have been obtained. To search for vacuum solutions of them in two spe- 
cific models, the conformal model and the zero-energy (Deser-Tekin) model, the dynamical systems have 
been derived, some de Sitter critical points and their stability have been investigated. These points are 
always exact constant solutions in the context of autonomous dynamical systems and describe the asymp- 
totic behavior. Some stable de Sitter critical points have been found. For any physical theories, to find 
exact mathematical solutions is an important topic. Next comes the physical interpretations of the solution 
thus obtained. Mathematically, de Sitter as the maximally space is undoubtedly important for any gravity 
theories. From observational side, recent studies illuminate that both the early universe (inflation) and the 
late-time universe (cosmic acceleration) can be regarded as fluctuations on a de Sitter background. So de 
Sitter takes a pivotal status in gravity, especially in modem cosmology. 

The solutions in section TV indicate that when f - 0, h ^ 0, the cosmological equations have stable de 
Sitter critical points. This means that the scalar ingredient h of torsion could be considered as a "phantom" 




I V I67 - 1 



(105) 
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field, since it does not interact directly with matter; it only interacts indirectly via gravitation. In the case 
f 0, h = 0, it is the pseudotrace axial ingredient / of torsion that produces the eff'ect of acceleration of 
cosmological expansion. Therefore the spacetime in the vacuum has the structure of de Sitter spacetime 
with torsion including the pseudotrace axial ingredient / as well as the scalar ingredient h. 

In summary, in the framework of gauge theory of gravity some cosmological models can be constructed 
to explain observable acceleration of cosmological expansion. The effect of acceleration of cosmological 
expansion in these models has the geometrical nature and is connected with geometrical structure of physical 
spacetime. The spacetime in the vacuum has the structure of de Sitter spacetime with torsion. 
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